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Abstract 

It is a review paper related to the following topics. General relativity (GR) is presented in the field 
theoretical form, where gravitational field (metric perturbations) together with other physical fields are 

>n : 

propagated in an auxiliary either curved, or flat background spacetime. A such reformulation of GR is 
exact (without approximations), is equivalent to GR in the standard geometrical description, is actively 
used for study of theoretical problems, and is directed to applications in cosmology and relativistic 
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astrophysics. On the basis of a symmetrical (with respect to a background metric) energy-momentum 
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tensor for all the fields, including gravitational one, conserved currents are constructed. Then they are 
expressed through divergences of antisymmetrical tensor densities (superpotentials) . This form permits 

(N ■ 

■ to connect a necessity to consider local properties of perturbations, which are analyzed in application 

oo 
in 

gauge invariance is studied, and its properties allow to consider the problem of non-localization in exact 

in ■ 

mathematical expressions. The M/string considerations point out to possible modification of GR, for 
example, by adding "massive terms" including masses of spin-2 and spin-0 gravitons. A such original 
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modification on the basis of the field formulation of GR is given by Babak and Grishchuk, and we 
present and discuss it here. They have shown that all the local weak-field predictions of the massive 
theory are in agreement with experimental data. Otherwise, the exact non-linear equations of the new 
theory eliminate the black hole event horizons and replace a permanent power-law expansion of the 
homogeneous isotropic universe with an oscillator behaviour. One variant of the massive theory allows 



j—j ' "an accelerated expansion" of the universe. 

1 Introduction 

Frequently in general relativity (GR) investigations are carried out under assumption that perturbations 
of physical fields propagate in a background given (fixed) spacetime, flat or curved, which is a solution 
to the Einstein equations [1] - [3]. A majority of tasks in modern cosmology and astrophysics are just 
formulated as the study of generation, propagation, evolution and interaction on different backgrounds. Exact 
cosmological and black hole solutions, for example, are used as such backgrounds. But, linear approximation 
(without taking into account "back reaction" ) , and flat or strongly simplified backgrounds frequently only 
are considered. Additional assumptions are used, but it is not clear how results will be changed without 
these assumptions, etc. All of these creates a necessity in a generalized (united) description of perturbed 
systems in GR. It is more natural and desirable to to present them in the framework of a field theory in a 
fixed background spacetime with the following properties: 
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• Constructions are covariant and give a possibility to use an arbitrary curved background — a solution 
to GR. 

• A perturbed system is defined by a) Lagrangian (action); b) perturbed (field) equations; c) conserved 
quantities, like energy, its density, etc. 

• Gauge (inner) transformations are defined explicitly with well described properties convenient in ap- 
plications. 

• There are no restrictions in approximations, thus perturbed equations, gauge transformations and 
conservation laws are to be exact. 

Both in theoretical analysis of perturbations and in applications, definitions of energy, momentum, an- 
gular momentum, their densities, and conservation laws for them turn out especially important. It is well 
known that conserved quantities, like energy, are not localized in GR. This means that it is impossible to 
construct unique covariant densities of these quantities in general. A reason is in physical foundations of 
GR, namely in the equivalence principle (see, e.g., the textbook [2]). In mathematical terms the problem 
is explained by a double role of a spacctimc in GR. On the one hand, it is an arena, on which physical 
fields interact and propagate; on the other hand, the spacetime itself is a dynamical object. Because the 
reason is objective, sometimes an opinion arises that the study of the conserved quantities in GR is senseless. 
The author of the paper, the same as many other researches, do not support this point of view. The fact 
of the non-localization cannot suppress the notions, like energy, themselves. Without doubts, gravitational 
interaction gives a contribution into a total energy, momentum, etc. of gravitating systems [2]. For exam- 
ple, to describe a binary star system one has to include a notion of gravitational energy as an energy of a 
connection. Considering gravitational waves in a bounded domain of empty space one can show that this 
domain has a total positive energy, etc. Such examples show that conserved quantities in GR are defined in 
non-contradictive manner, at least, as non-local characteristics, and, of course, have to be examined. 

Total energy-momentum, angular momentum of an asymptotically flat spacetime have been studied in 
details and are continuing to be considered. One of the main successes was the proof of the positive energy 
theorem [4]. This initiated a renewed interest to the problem. Energy-momentum and angular momentum 
became to be associated with finite spacetime domains. Such quantities are called as quasi-local ones and their 
examination during last two-three decades was very successful (see a recent nice review [5] and references 
there in). Returning to the cosmological problems, where conversely more frequently local properties of 
perturbations are studied, we accent a necessity to connect local and quasi- local derivations. 

There are many of approaches in GR where both an evolution of perturbations and conservation laws 
for them are studied. In this paper we consider only one of approaches, which in a more measure satisfies 
the above requirements. The perturbed Einstein equations are rewritten as follows. The linear in metric 
perturbation terms are placed on the left hand side; whereas all nonlinear terms are transported to the 
right hand side, and together with a matter energy-momentum tensor are treated as a total (effective) 
energy-momentum tensor t^"^. This picture was developed in a form of a theory of a tensor field with 
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self-interaction in a fixed background spacetime, where t^ is obtained by variation of an action with 
respect to a background metric. Frequently it is called as a field theoretical formulation [6] of GR, we 
will call it simply as the field formulation, and we just review it here. The active history of these studies 
was begun in 50-th of XX century. Deser [7] has generalized previous works and suggested the more clear 
presentation without expansions and approximations in Minkowski spacetime. We [8] developed the field 
formulation of GR on arbitrary curved backgrounds. Advantages of a such description was demonstrated 
in several applications. A closed Friedmann world was presented as a gravitational field configuration in 
Minkowski space [9]; trajectories of test particles at neighborhoods of event horizons of black holes were 
analyzed [10]; t^^f 1 and its gauge properties were used for a development of elements of quantum mechanics 
with gravitational self- interaction [11], in the frame of which some of variants of an inflation scenario were 
examined [12]; a distribution of energy in black holes solutions was constructed [13]; it was studied so-called 
the weakest falloff conditions for asymptotically flat spacetime at spatial infinity [14]. A related bibliography 
of earlier works particularly can be found in [6] - [8]. The foundation of the field formulation of GR and 
some references can be also found in the review and discussion works [15, 16]. A more full, as we know, the 
modern bibliography can be found in the works [17]. 

There are different possibilities to arrive at the field formulation of GR. Deser [7] used a requirement: • 
for a linear massless field of spin 2 in a background spacetime it has to be a source in the form of a total 
symmetrical energy-momentum tensor of all the fields including gravitational one. Namely this principal 
was used as a basis for constructions in [8]. It is well known the other method, which was more clearly 
represented by Grishchuk [6] and which briefly is formulated as a transition • from gravistatics (Newtonian 
law) to gravidynamics (Einstein's equations). Keeping in mind gauge properties of the Einstein theory 
one obtains the field formulation of GR as a result of • a "localization" of Killing vectors of a background 
spacetime, which have a sense of parameters in an action of a gauge theory [18]. The way, which has an 
explicit connection with the standard geometrical formulation of GR is based on a simple • decomposition 
of usual variables of the Einstein theory onto dynamical and background quantities [19]. 

In section 2, a construction of the field formulation is given on the basis of the last of the above methods. 
In section 3, using the results of the works [20, 21] we present conservation laws in the field formulation of GR. 
The conserved currents are constructed on the basis of symmetrical energy-momentum tensor and express, 
thus, local characteristics of conserved quantities. At the same time the currents are derived as divergences 
of antisymmetrical tensor densities (superpotcntials), integration of which just leads to surface integrals, 
which are quasi-local conserved quantities. In section 4, we give a generalization of the results of sections 
2 and 3 for various definitions of metrical perturbations and resolve related ambiguities. The one of more 
desirable properties is that an energy-momentum complex of a theory has to be free of the second (highest) 
derivatives of the field variables. The energy-momentum tensor in [8] does not satisfy this requirement. 
Babak and Grishchuk recently improved this situation [22] . Developing the approach [8] they reformulated 
the field interpretation of GR satisfying the above property; in section 5 we outline their approach and 
details of the method. Gauge transformation properties in the field formulation of GR and their connection 
with the non- localization problem in GR are discussed in section 6. The original and perspective technique 
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developed in work [22] is naturally generalized for constructing a gravitational theory with gravitons of non- 
zero masses in the work [23]. In its framework Babak and Grishchuk have also found and examined static 
spherically symmetric solutions in vacuum, as well as homogeneous and isotropic cosmological solutions. All 
of these are included in section 7. We discuss some questions in the last section 8. Thus, the goal of the 
present review is a presentation of the field approach to GR with its development permitting to arrive at the 
Babak-Grishchuk massive gravity. This gives a possibility to understand better the properties of this version 
of the massive gravity and its connection with GR. Here, we do not consider other variants of theories with 
massive gravitons. 

Below we give more general notations used in the paper. 

• Greek indexes mean 4-dimensional spacetime coordinates. Small Latin indexes from the middle of 
alphabet i, j, k, . . ., as a rule, mean 3-dimcnsional space coordinates; large Latin indexes A, B, C, ... 
are used as generalized ones for an arbitrary set of tensor densities, like Q A . Usually x° — ct, where 
c is speed of light; k — 8nG/c 2 is the Einstein constant; {af}) and [a(3] are symmetrization and 
antisymmetrization in a and /3. 

• The dynamic metric in the Einstein theory, as usual, is denoted by g^ u (g — detg^ u ), whereas ~g 
(7j = det g^) is the background metric; t]^ is a Minkowskian metric. A hat means that a quantity "Q" 
is a density of the weight +1, it can be Q = \J~^gQ, or Q = \/~^§Q 1 or independently from metric's 
determinants, it will be clear from a context. A bar means that a quantity "Q" is a background one. 
Particular derivatives are denoted by {pi) and (d a ); (D a ) is a covariant derivative with respect to g^ v 
with the Chistoffel symbols T^ 7 ; (D a ) is a background covariant derivative with respect to ~g^ v with 
the Chirstoffel symbols T^ 7 ; 6/SQ A — the Lagrangian derivative; £^Q A = -£ a D a Q A + Q A \Z D a ^ 

— the Lie derivative of a generalized tensor density Q A with respect to the vector ; for example for 
the contravariant vector one has £^Q^ — —£, a D a Qv + D a ^. 

• R a nf3v, R^v, G>„, and R are the Riemannian, Ricci, Einstein and matter energy-momentum 
tensors and the curvature scalar for the physical (effective) spacetime; R° i? M „, G^, T^ u and R 
are the Riemannian, Ricci, Einstein and matter energy-momentum tensors and the curvature scalar 
for the background spacetime. 

2 Exact perturbed Einstein equations on an arbitrary curved back- 
ground 

At first we define the Lagrangian for the perturbed system. Consider the usual action of GR: 

S=-J d 4 xC E = - JL J d 4 xR(g^) + -J d 4 x£ M ($ A , g^) (2.1) 
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where for the sake of simplicity one assumes that £ M (<J> A , g^) depends on the first derivatives only. Let us 
write out the Einstein equations together with the matter ones in the form: 

SC E 1 SR 6C M n 

— — = — — — T7 h — — = , 2.2 

bg^ v 2n5g^ Sg^ u v ' 

Now, define the metric and matter perturbations as 

V=gg'"' = g^ =W + P V , <& A = <H + 4> A . (2.4) 
The background system is described by the action: 

S = - c J d^xjE = __L J d i x R + ij £ x Jm i (2.5) 

and the background quantities g^ v and & A satisfy the corresponding background equations: 

1 SR Si™ si™ 

+ T^ = ' T^ = °- (2-6) 



2k 5g^ v Sg^ v S$ A 

The perturbations and <ft A now are thought as independent dynamic variables. The perturbed system 

is to be described by a corresponding Lagrangian on the background of the system (2.5) and (2.6). Let us 

construct it. Substitute the decompositions (2.4) into the Lagrangian of the action (2.1), subtract zero's and 

linear in and (f> A terms of the functional expansion, and add any divergence: 

xrE xrE i i 

C dyn = c E {g + l, $ + 4>) - l^-^= - <p A -^= - C E - —d a k a = ~—C 9 + t m , (2.7) 

Sg^ v 8§> A 2k 2k 

see on functional expansions in the the book [3]. Lagrangians, like (2.7), are called as dynamical Lagrangians 

in the terminology of [19]. Zero's term is the background Lagrangian, whereas the linear term is proportional 

to the l.h.s. of the background equations (2.6). However, one should not to use the background equations in 

£d V n b e f ore it s variation because really C dyn is not less than quadratic in the fields and cf) A in functional 

expansions. 



If one chooses the vector density 



with the definition 



k a ee g""A^ - g^A^ (2.8) 



A°„ = - = \g a » (D^g pu + D v g pp . - D p g^) , (2.9) 
where the decomposition (2.4) is used, then a pure gravitational part in the Lagrangian (2.7) is 



C g = R{g^ + l^)-l^R^-g^R^ + d^ 

= -(A£„ - A°JZ)D p l^ + {W + in (A^A^ - A^A^) . (2.10) 

It depends on only the first derivatives of the gravitational variables . In the case of a flat background 
the Lagrangian (2.10) transfers to the Rosen covariant Lagrangian [24]. The matter part of the dynamical 
Lagrangian (2.7) is 

- Sf M Sf M — — 
C m = C M (g + l , * + <f) -r=-/^=-£«. (2.11) 
V ; 5g^ 53> A 
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The variation of an action with the Lagrangian C dyn with respect to and some algebraic calculations 
give the field equations in the form: 

+ #£„ = « (t»„ + C) = . (2-12) 

where the l.h.s. linear in and </> A consists of the pure gravitational and matter parts: 

= ^l pa -= = \ (D P D% V + g^D p Dj»° D p DjP D P D P IJ>) , (2.13) 



The r.h.s. of Eq. (2.12) is the total symmetrical energy-momentum tensor density 

xrdyn £ / I \ 

4' ot) = 2-=- = 2—= -— £ 9 + L m = i» + . (2.15) 

In expansions, is not less than quadratic in f^ 1 " and yl that follows from the form of the Lagrangian 
C dyn . The explicit form of the gravitational part is 

1 



t 9 

' xv ~ K 



Sffi + l9^9 pa ) (A^A^ - A^A^) + D T QlJ\ , (2.16) 



where 



2Q^ = -9^i af) Kp + huK ( fg af) - - il^a + i 0T (^IpSau + K f3 g ap ) 

+ K (K, ~ A^«T<w) + # (A^ - A^g a ,) . (2.17) 

The matter part is expressed through the usual matter energy- momentum tensor of the Einstein theory 

as 

= [W - l9^9 pa ) {T P a - h9 P *T„x9 nX ) - T pv \ 

-2 * ff^+^^j, (2.18) 

and is also not less than quadratic in and 4> A in expansions. At the usual description of GR the definition 
of the energy-momentum tensor by S£ E /6g^ v is senseless because it is vanishing on the Eq. (2.2), whereas 
i^"^ defined in (2.15) is not vanishing on the field equations (2.12). A formal reason is that in the Lagrangian 
(2.7) the linear terms are subtracted. 

By the definitions (2.14) and (2.18), the field equations (2.12) can be rewritten in the form: 

&j; v = K(% v + 6i%,)=K%p (2.i9) 

where Siff^ = — iff„ is equal to i™ v in Eq. (2.18) without the second line. Thus 8t^f v can be thought as 
a perturbation of i*f v = T pui which includes even linear perturbations in the dynamic fields and does not 
follow from the Lagrangian (2.7). But now ijf/^ is the source of the linear gravitational field only without 
linear matter part (see Introduction). 
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Let us demonstrate the equivalence with the Einstein theory. Transfer i^"^ to the l.h.s. of Eq. (2.12) 
and use the definitions (2.13), (2.14) and (2.15) with (2.7): 



r; L -l<t> l — K f 



(tot) 



-2k 



dgP° 5 



+ 2k 



S 



~R (g al3 + Z Q/3 ) + C M ($ A + (f> A ; g»» + / >y ) 



-—R + C M } . 

2k 



(2.20) 



Because the third line is proportional to the operator of the background equations in (2.6), then Eq. (2.12), 
as seen, is the Einstein equations (2.2), only in the form with using the decompositions (2.4). 



3 Conservation laws 

At the beginning we discuss differential conservation laws on Ricci-flat (including flat) backgrounds. One has 
to take into account Q A = 0, C M = 0, = and use SR/5g fJ - 1 ' = instead of the background equations 
(2.6). Then the Lagrangian (2.7) is simplified to 

C dyn = -^-C 9 + t m = -^-C 9 + t M (V; W + l^), (3.1) 
2k 2k \ J 

and the field equations (2.12) transform into the form of Eqs. (2.20): 

&^ = K{% v + fy,)=Kfft? ) - (3-2) 

Thus, for Ricci-flat backgrounds and i^J^ coincide, and Eqs. (2.20) and (3.2) have the form announced 
in Introduction. Because in Eq. (3.2) D„G^ U = 0, a divergence of Eq. (3.2) leads to 

Dj^ ot > = . (3.3) 

A contraction of with background Killing vectors A Q gives a current J V {X) = i$ ot ^ v \ li , which is also 
differentially conserved: 

D V J V (A) = %J V (A) = . (3.4) 

Integration of J v leads to non-local conserved quantities, Consider a background 4-dimensional volume 
Vi, the boundary of which consists of timelike "surrounding wall" S and two spacelike sections: So : = to = 
const and Si := t\ = const. Because the conservation law (3.4) is presented by the scalar density it can be 
integrated through the 4-volume V4: J v d^jP(X)d 4 x = . The generalized Gauss theorem gives 

f J°(\)d 3 x- f J°(X)d 3 x+ I J»(\)dS„ = (3.5) 

J Si J So J S 

where dS^ is the element of integration on S. If in Eq. (3.5) § s JP{\)dS tl = , then the quantity 

V{\) = J J°(\)d 3 x (3.6) 

is conserved on £ restricted by <9£, intersection with S. In the converse case, Eq. (3.5) describes changing 
the quantity (3.6), that is its flux through <9S. It can be also assumed dT, — > 00. 
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The differential conservation laws (3.3) have also a place for backgrounds presented by Einstein spaces in 
Petrov's defenition [25]: R pu = Ag where A is a constant (see [8, 26, 27]). For arbitrary curved backgrounds 
there are no conservation laws, like (3.3). Indeed, in the general case D v (g^ v + ^ in (2.12), and 

D V G^ 7^ in (2.19). The reason is that the system (2.7) interacts with a complicated background geometry 
determined by the background matter fields & A . Many of cosmological solutions are just not of the Einstein 
spaces. 

Conservation laws for arbitrary curved backgrounds and arbitrary displacement vectors £ Q were con- 
structed in [21]. With using the technique of canonical Noether procedure developed in [28] and applied to 
the Lagrangian (2.7) it was obtained the identity: 

-G^C + -l^RxvC + ( p = Dj» v = dj pu . (3.7) 

K K 

The superpotential has the form: 

P v = -F^D P C ] + V^x^ = - (l p[p D p C ] + &D a t^ - D^tX a ) (3.8) 

and, thus, d pv I pv = 0. It generalizes the Papapctrou superpotential [29]; indeed for the translations in 
Minkowski space £ A = #K in the Lorentzian coordinates one gets 

i(p) = = Y K d ° ( Spta - s P itMT - + r " ip ) • (3 - 9) 

The same superpotential (3.8) was constructed in [20] by the other way, namely, by the Bclinfantc sym- 
metrization of the canonical system in [28]. The last term on the l.h.s. of (3.7) is 

2kC = 2 (z pa D p l p - i pa D p z%) - (z^D^ - i pa D^z p ^ + {l pv D v z - zDJ p ^ (3.10) 

where 2z pa = —£^g pa , and, thus, disappears on the Killing vectors of the background. 

To write out a physically sensible conservation laws from the identity (3.7) one has to use the field 
equations, which we substitute in the form of Eq. (2.19). Then the identity (3.7) transforms to the equation 

I" = TfC + C M = Dj^ = dJ plJ . (3.11) 

The generalized total energy-momentum tensor density is 

f p = tl p + St?" + -l pX R Xu = + -l pX R Xu (3.12) 

K K 

where i p tot ^ v is exchanging with , and the interaction with the background geometry term, l pX R\ l/ , is 

adding. Thus, T p plays the same role as i^ ot ^ v in Eq. (3.4) if Killing vectors exist. However, the current 
I p has a more general applicability, than J p in (3.4): it is conserved, D p I p = d^I p = 0, on arbitrary 
backgrounds and for arbitrary £ Q . It is important, for example, for models with cosmological backgrounds 
where not only the Killing vectors are used fruitfully (see, e.g., [30]). 

Due to antisymmetry of the superpotential (3.9) the conserved quantity, like (3.6), is expressed over a 
surface integral 

V(0 = I i 0k (0ds k (3.13) 

JdY. 
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where dsk is the element of integration on 9S. It is important expression because it connects a quantity 
■P(£) obtained by integration of local densities with a surface integral playing a role of ^quasi-local quantity 
(see discussion in Introduction). 



4 Different definitions for perturbations 

In GR, components of each of metrical densities 

9 a = {<T, 9^ y/=9<r, V=99^, (-ff)<T, ■ ■ •} (4.1) 

could be chosen as independent dynamic variables. In the terms of generalized variables (4.1) the action of 
GR (2.1) is rewritten as 

S = - [ al A xt E ^ = — f d 4 xR(g a ) + - [ d 4 xl M ($ A , g a ) . (4.2) 

c J 2kc J c J 

Variation with respect to g a gives the gravitational equations in a corresponding form instead of (2.2). The 
perturbations could be also defined for each of metric variables in Eq. (4.1) as 

{ 9 a =r + h a } = [ 9ia> = g^ + v, gr = ¥ s + i"", <r = r v + n™, . . .} . (4.3) 

For the decomposition (4.3) following to the rules of constructing the Lagrangian (2.7) one gets 
%d v n — _J_p^_L(,^j_r M /'iI_L^. 



C7T = --R (g^ + h a ) + C M [$> A + <j> \ Y+h a 



- ^ + ^ \-^J=A-^-R + l^\-Uk». (4.4) 



2nSg a Sg a J S$ A \ 2n J 2k 

Its variation with respect to h a and some re-calculations give the Einstein equations in the form (2.12): 

QL{a) + ^L(a) = ftot a) u ^s 

The total symmetrical energy-momentum tensor density is defined as usual: 

ftct a) = 2 %= . (4.6) 
In the l.h.s. of Eq. (4.5) the independent variables h a are replaced by the other variables 

l (a)- n ffgS> ^ J > 

which are also considered as independent ones due to the background equations. Thus, the same operators 
(2.13) and (2.14) are applied to l 1 ?" 

For some of different decompositions (4.3): g\ — gf + hi and g% = g% + h% the variables (4.7) differ 
one from another in the second order in perturbations: 1^2) = ^(Ti) + P(a)\2- Because differences inter the 
linear expressions of equations (4.5) the energy-momentum tensor densities a1 ^ and a2 ^ get the same 
differences too. Firstly, for the case of flat backgrounds this fact was noted by Boulware and Deser [31]. 
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For the system (4.4) the identity 

takes a place and has exactly the form of the identity (3.7) with replacing with Z^ only. Substituting 
Eq. (4.5) into the identity (4.8) we obtain the conservation law in the form: 

iU - + < (a) " + T + C? o) = %X + Cf o) = (4.9) 

analogous to (3.11) with (3.12). Thus a family of conservation laws (4.9) presents a corresponding family of 
superpotentials, which can be presented in the form of the Abbott-Deser type [26]: 

% = \ (l&D^ + &DJZ) D^l^) . (4.10) 

It is exact the formula (3.8) with exchanging by Indeed, the known Abbott-Deser superpotential 
inter just this family with the decomposition g pv = g + h^ v through the transformation (4.7) and for 
(anti-)de Sitter's backgrounds. Superpotentials of this family differs one from another due to the difference 
in perturbations. Otherwise, there is an independent way. Thus in the work [20], the generalized Belinfane's 
technique has been applied for "symmetrization" of the conserved quantities suggested by Katz, Bichak and 
Lynden-Bcll [28]. Their method does not depend on the choice of the variables from the set (4.1) and gives 
uniquely [21] the superpotential (3.8). Thus, theoretically on the level of superpotentials the Boulware-Deser 
ambiguity [31] is resolved in favor of l pv . 



5 Gravitational energy-momentum tensor 

In this section we present the results by Babak and Grishchuk [22]. Following them here and in next 
sections we consider the Minkowski space with R ap f3a- — as the background spacetime. Our presentation 
is technically simpler than that of Babak and Grishchuk [22], although, of course, is equivalent to their 
one. Thus we repeat their calculations on the basis of formulae (2.7) - (2.20) simplified to the case of Eqs. 
(3.1) - (3.3). Also, in the work [22] as independent variables it is used = l pv j \/—g, whereas we use 
V 1V . Leaving the field equations to be equivalent, this leads to a difference in the direct definitions of the 
energy-momentum tensors. However, with taking into account the field equations this difference disappears 
and does not influence on results and conclusions. 

With using the definition (2.4) we present the expression (2.9) through the gravitational variables l^ u : 



A* = -L- 



9 PP DJ xp + g vp Dj x P p9 Xp D p l 



a/3 



+ | (g aP 8 x Dj^ + g afj 5 x v Dj^ - g a0 g^g xp D p i a ^] (5.1) 



where g pu , g^ u and \J~-g arc thought as dependent on the definition (2.4). Substituting Eq. (5.1) into Eq. 
(2.16) with (2.17) one finds that depends on the second derivatives of l pv ' . After using the field equations 
(3.2) the second derivatives are left, but only minimally, as 



C = C-d) + Q af3 ^(Cp ~ + VyFlr^D^W* - W); (5.2) 

(V^v) 2 Q afhll/ = i a(lt T )f> + i p{M f )a + i a(il i v)0 - Wi af> - (Jf 3 + i af> ) ■ (5.3) 
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The reduced part with only the first derivatives is 

1 



•//.// 



2D p l^Dj pa - 2Dj^ a D p r> 3 



(g-red) 4K^f 

- lg 0p g a ^Dj^Djr° 

+ l(2g^g^-g^g u5 )(2g pa g^-g a pg pa )D 5 i pa Dj^] . (5.4) 

The matter part in (5.2) has appeared due to using the field equations (3.2). 

In [22] it was suggested the original way to exclude the second derivatives from the energy-momentum 
tensor without changing the field equations. The Lagrangian (2.10) was modified as follows 

£ 9 (mod) =£ 9 + ^ Ppa R apf 3<r. (5.5) 

This is a typical way of incorporating constraints (because R ap j} a = 0) by means of the undetermined 
Lagrange multipliers. The multipliers A al3pa form a tensor which depends on ~g pv and V" 1 (without their 
derivatives) and satisfy A a/3pCT = -AP siaa = -A"' 7 '" 3 = A^p. Thus, the field equations (3.2) do not 
change. Then, in a correspondence with the modified Lagrangian (5.5), the modified energy-momentum 
tensor density is 

d \Zd) = k£9 - D af 3 (A^ Q/3 + A 1 ^) (5.6) 

instead of (2.16). The originally undetermined multipliers A p " a P will now be determined. They can be chosen 
in a such way that the remaning second derivatives in (5.2) can now be removed. The unique possibility is 
AP ua P = {l av pP - /4^/=f . Thus the equations (3.2) arc not changed, but they have to be rewritten 

in the form 



a/3 



(g^ + l» u )(g a P + l afi ) - (gP a + l» a )(g»P + l ufi ) 



- K (*g(mod) + *m ) ~ ^(mod-tot)' ( 5 ' 7 ) 

Here, the r.h.s. defined as a symmetrical (metric) energy-momentum tensor for the system (5.5) is the source 
for the generalazed d'Alembert operator (general wave operator). Thus the l.h.s. in (5.7) is not more linear 
in l pv . Because on the flat background the divergence of the l.h.s. in (5.7) is identically equal to zero, 
then D u i^ iod _ tot ^ = 0. In Eqs. (5.7), i^ mod -tot) can be reduced by the equations of motion, then they are 
rewritten as 



L{mod) 



1 (g-red) + ^ \ l afi 2$ afi 1 p I T z y,v ~ Kl (mod-tot-red) \°-°) 



Thus, indeed on the equations of motion the energy-momentum tensor density in (5.8) is only with the first 
derivatives of gravitational variables. Again D lJ t^ od _ tot _ red ^ = 0. Let us show that Eq. (5.8) is equivalent 
to the usual Einstein equations. Multiplying it by \/—g, and using the identification (2.4), the definition 
(2.18) for the flat background and the definition (5.3), in the Lorentzian coordinates, one easily gets 

\d a0 [(-g)(g^g a0 - g» a g^)] - k(- 9 ) (t? L + T^) . 
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After substituting the Einstein equations nT^ v = C" this equation transfers to the identity. Thus, indeed 
Eq. (5.8) is equivalent to the Einstein equations, and one finds that {—g)t^ L is the Landau-Lifshitz's 
pseudotensor [1]. After all, one concludes that tf" ed \ is the covariantized pseudotensor (— g^lW / V~V ■ 



6 Gauge invariance properties 

Properties of the field formulation of GR under gauge transformations follow from the usual covariant 
invariance properties of GR in the geometrical description. We demonstrate it briefly (for details see [8, 19]). 
Consider the same solution to GR ^^(x) and g' lll '{x l ) presented in two different coordinate systems: {x} 
and {x'} connected by the coordinate transformation x' = x'(x). Now let us do a decomposition of the 
type (2.4) in both the cases: g^ v {x) = g~^{x) + l^{x) and g ,fl "(x') = g^{x') + i'> iV {x'), with the same 
form of the background metric gv. For the solution in primed coordinated from the points with coordinate 
quantities x' one displaces to the points with coordinate quantities x. After that one has to compare both 
solutions. An analogous procedure has to be carried out under the matter variables, assuming the coordinate 
transformation in the form: 

x , a = xa + r + + ±_ e (f^) p + . . . , 

where £ M are assumed as enough smooth, we obtain the transformation [8, 19]: 

= j> + - 4 + , = ^ + _ £ \ ($ A + ^) , (6.i) 

k=l k=l 

which are called as gauge (inner) transformations in the field formulation of GR. Indeed, they do not affect 
both the coordinates and the background quantities. 

Firstly let us consider the field formulation of the section 2. It is not difficult to see that the Lagrangian 
(2.7) is invariant under the transformation (6.1) up to a divergence on the background equations (2.6). The 
expression (2.20) gives a possibility to understand that equations (2.12) are gauge invariant on themselves 
and on the mentioned background equations. But the energy-momentum tensor density (2.15) (or (2.19)) is 
not gauge invariant. Even on the field equations one has 

< ( t ot) = 4^ + ^ (*' - o + {I'-utf-®, 

>Jfc m = Kiff»+G^(l'-l). (6.2) 

A 

The transformations (6.1) with $ = 0, as it has to be for i? M „ Q/ 3 — 0, are also gauge transformations for 
the formulation in section 5. The Babak-Grishchuk Lagrangian is also gauge invariant up to a divergence, 
and the equations (5.8) are gauge invariant on themselves. Concerning the energy-momentum tensor, on the 
field equations one has 

^(mod-tot-red) = ^(rnod-tot-red) + QL(mod) _ ^ (g 3) 

The non-localization problem of energy and other quantities in GR became evident from the moment 
of constructing GR beginning from the original Einstein's works. During prolonged time it was illustrated 
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by a non-covariance of pseudotensors. The use of an auxiliary background spacetime permitted to consider 
covariant conserved quantities, but the non-localization became to be explained by an ambiguity in a choice 
of a background. However, on this level there was no suggested an unique mechanism for description of this 
ambiguity. The use of the gauge transformation properties in the field formulation of GR closes this gap in 
GR. At the beginning of this section we just express the connection between different choices of backgrounds 
explicitly. A gauge non-invariance in the energy-momentum tensors (6.2) and (6.3) just expresses the non- 
localization of energy, momentum, etc. in GR in exact (without approximations) and explicit mathematical 
expressions. It is a one of advantages in using the field formulation of GR. 



7 Gravity with non-zero masses of gravitons 

Babak and Grishchuk using their technique [22] have constructed a variant of theory of gravity with non-zero 
masses of gravitons with interesting properties [23]. Following to [23] independent variables V 1 " — ' j 'y/—g 
are used. It is natural to assume that the Lagrangian may also include an additional term similar to the one 
in Eq. (5.5), but where the quantity R ap fja is the curvature tensor of an abstract spacetime with a constant 
non-zero curvature: R ap (3a = K (jjupgpa — 9aa9p(3) where K is with the dimensionality of [length] -2 . If one 
adds k af) P a R ap[ja with A^ a/3 = {A^g)- 1 (i av ft» - Hz>^ , changing g» v -> g^ , then the additional term 
in the Lagrangian (5.5) is \\J— gK {l a ^l a (3 — ^^a) • Clearly, the new theory is not GR, but one recognizes 
in this term the Fierz-Pauli mass-term [32]. Thus, noting that the structure (5.5) generates mass terms 
and finding that only two independent quadratic combinations of V 111 exist, Babak and Gishchuk arrive at a 
2-parametric family of theories with the additional mass terms in the gravitational Lagrangian (5.5): 

^mass) = 4^) + [kiltie? + W] , (7.1) 

where k\ and ki have a dimensionality of [length] -2 . 

Of course, the additional term in (7.1) gives a contribution both into the r.h.s. and into the l.h.s. of Eq. 
(5.8), and the equations of the new gravity theory symbolically could be rewritten as 

= K tf» . . (7.2) 

L(mass) {tot— mass) v ' 

These equations are, of course, covariant, however, unlike Eq. (3.2) and Eq. (5.8), the new field equations 
(7.2) are not gauge invariant. There arc no transformations, like (6.2). Thus, there is no a problem with a 
localization of vfY t n — it is localized! 

(tot— mass) 

To have a direct comparison with GR effects it is more convenient to present Eq. (7.2) in the quite 
equivalent quasi-geometrical form: 

G pv + M„„ = kT„„ (7.3) 

where the massive term is 

M„„ = {25°&Z - g a!i g pv ) (hl a0 + k 2 g af) F p ). 

Notice that the Bianchi identity D U G^ = takes a place in effective spacetime. Besides, with taking into 
account the matter equations (2.3) one has D„T^ — 0, as usual. Thus, after differentiation of Eq. (7.3) 
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one obtains D u M plJ = 0. Although these equations are merely the consequences of the full system (7.3), 
and therefore contain no new information, it proves convenient to use them instead of some members of the 
original set (7.3). 

To give a physical interpretation of k\ and k 2 , following to the analysis by Ogievetsky and Polubarinov 
[33], and by van Dam and Veltman [34], one considers the linearization of the Eqs. (7.3): 

\ (D p D P l^ + g plJ D p DJ p ° - DpDJP - D p D p lf) + 2k Y l pv - (h + 2k 2 )g p J a a = 0- (7.4) 

The divergence of this equation is 

D v [2hF v - (h + 2fc 2 )^X] - 0, (7.5) 

which is the linearized version of the equation D v M pv = 0. 

Consider the first case with k\ ^ k 2 . The full system (7.4) is equivalent to 

aH^ + a 2 H" u = 0, (7.6) 
Ul a a +pH a a = 0, (7.7) 

together with Eq. (7.5). Here, □ = g afll D a f3, 



ttv-v = h ^ _ h+k 2 _ k x +k 2 — P u fci +k 2 



with g^ v H^ v = and D V H^ V = 0. Thus, parameters in the wave-like equations (7.6) and (7.7) are 



a^4 kl , 2 = _ 2k 1 (k 1+ Ak 2 ) ^ 

h + k 2 

They can be thought as inverse Compton wavelengths of the spin-2 graviton with the mass m 2 — aU/c 
associated with the field and of spin-0 graviton with mass m = (3h/c associated with the field 

With studying the weak gravitational waves in the massive gravity one finds certain modifications of 
GR. Thus the spin-0 gravitational waves, presented by the trace I™ = l a,3 r] a j3, and the polarization state of 
the spin-2 graviton presented by the spatial trace H lk r]ik both, unlike GR, become essential. They provide 
additional contributions to the energy-momentum flux carried by the gravitational wave, and the extra 
components of motion of the test particles. However, gravitational wave solutions, their energy-momentum 
characteristics, and observational predictions of GR are fully recovered in the massless limit a — > 0, (3 — > 0. 

For the case with the mass term of Fierz-Pauli type, fci + k 2 = 0, that corresponds 1 — > oo (see (7.9)), 
the full set of equations (7.4) is equivalent to 

l a a = 0, Ol^ + 4fciP y = 0, D v l» v = 0. 

This case is interpreted as unacceptable [23]. Even in the limit a — > 0, there remains a nonvanishing 
"comoving mode" motion of test particles in the plane of the wave front. The extra component of motion 
is accounted for the corresponding additional flux of energy from the source, typically, of the same order of 
magnitude as the GR flux. This, at least, is in a conflict with already available inderect gravitational- wave 
observations of binary pulsars [35]. Such theories probably have to be rejected. 
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In [23], the full non- linear equations (7.3) were analyzed from the point view of the black hole and the 
cosmological solutions. Thus, searching for static spherically-symmetric solutions in vacuum it is necessary to 
consider three independent equations from (7.3), unlike GR where there are two ones only. The consideration 
is simplified if one assumes a = f3, however all the qualitative conclusions remain valid for a ^ (3. Combining 
analytical and numerical technique Babak and Grishchuk have demonstrated that the solution of the massive 
theory is practically indistinguishable from that of GR for all 2M <C R <C 1/a, where R and M are the 
radial and mass parameters of the Schwarzschild solution. For R larger than 1/a the solution takes the form 
of the Yukawa-type potentials; therefore they call this massive theory as finite-range gravity. The solution 
of new theory is also deviate strongly from that of GR in the vicinity of R — 2M that is the location of 
the globally defined event horizon of the Schwarzschild black hole in GR. In the massive gravity the event 
horizon does not form at all, and the solution smoothly continues to the region R < 2M and terminates 
at R = where the curvature singularity develops. Since the aM can be extremely small, the redshift of 
the photon emitted at R — 2M can be extremely large, but it remains finite in contrast with GR solutions. 
Infinite redshift is reached only at the singularity R — 0. In the astrophysical sense, all conclusions that rely 
specifically on the existence of the black hole even horizon, are likely to be abandoned. It is very remarkable 
and surprising that the phenomena of black hole should be so unstable with respect to the inclusion of the 
tiny mass-terms, whose Compton wavelength can exceed, say, the present-day Hubble radius. 

It was also considered homogeneous isotropic solutions in the framework of the massive gravity. Matter 
sources were taken in the simplest form of a perfect fluid with a fixed equation of state. There are two 
independent field equations from the set (7.3), unlike GR where there is only one in the same case. First, 
if the mass of the spin-0 graviton is zero, 1 — 0, the cosmological solutions are exactly the same as those 
of GR, independently of the mass of the spin- 2 graviton, i.e., independently of the value of a 2 . This result 
is expected due to the highest spatial symmetry: the spin-2 degrees of freedom have no chance to reveal 
themselves. Then, for (3 2 ^ it was considered technically more simple case A0 2 — a 2 which was studied 
in full details. Qualitative results are valid for A0 2 ^ a 2 . Again, combining analytical approximations and 
numerical calculations it was demonstrated that the massive solution has a long interval of evolution where 
it is practically indistinguishable from the Friedmann solution of GR. The deviation from GR are dramatic 
at very early times and very late times. The unlimited expansion is being replaced by a regular maximum of 
the scale factor, whereas the singularity is being replaced by a regular minimum of the one. The smaller (3, 
the higher maximum and the deeper minimum, i.e., the arbitrary small term in the Lagrangian (7.1) gives 
rise to the oscillatory behaviour of the cosmological scale factor. 

Following the logic of interpretation that a 2 and (3 2 define the masses, they are thought as positive. 
However, the general structure of the Lagrangian (7.1) does not imply this. Then, if one allows a 2 and (3 2 
to be negative, the late time evolution of the scale factor presents an "accelerated expansion" that is similar 
to the one governed by a positive cosmological A-term. The development of this point could be useful in the 
light of the modern cosmological observational data [36] . 
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8 Concluding remarks 



Due to the non-localization problem/property a definition and a study of conserved quantities in GR are 
not trivial tasks. Then rather one has not to follow unconditionally to some one unique method. However, 
to restrict such methods, one has to examine their possibilities to satisfy known natural tests. As a rule, in 
applications of expressions of such approaches it is required: a) the energy density for the weak gravitational 
waves has to be positive [1]; b) ratio of mass to angular momentum in the Kerr solution has to be standard 
[28]; c) one has to obtain the standard conserved quantities both at spatial and at null infinity for asymp- 
totically flat solutions [2]. The usual field formulation (UFF) of GR presented in sections 2 and 3, and the 
Babak-Grishchuk modification presented in section 5, they both satisfy all these evident requirements, and 
one cannot do a choice on this basis. 

In Introduction it was presented reasons why it is important to study perturbations on arbitrary curved 
backgrounds, and to present conserved currents as divergences of superpotcntials. In particular, just these 
requirements initiated a development of UFF. Conversely, the opinion of Babak and Grishchuk is that it is 
enough to use only the background Minkowski space. This has its own convincing foundation. All the modern 
"direct" experiments, as well as a aforementioned theoretical restrictions a) - c), use as a background a flat 
spacetime. Besides, the field formulation even can describe arbitrary curved and topologically non-trivial 
solutions of GR as a field configuration in Minkowski space (see [9, 10]). Paying a necessary opinion to the 
position of these authors, we note that there are no principal arguments against a following generalization 
of the modified Babak-Grishchuk formulation of GR onto arbitrary curved backgrounds and constructing 
conservation laws with the use of superpotentials. 

Let us return to the requirement to have only the first derivatives in the energy-momentum tensor, one of 
the main reasons of which is a correct formulation of the initial problem. However, at least, from this point 
of view second derivatives, which appear in the framework of UFF do not initiate a criticism. One needs 
to consider the energy-momentum tensor density i 9 and the conserved current 1^ (see (2.10) and (3.11) 
with (3.12)). In [21] we have shown that the symmetrized quantities constructed in [20] coincide with the 
ones presented here and included into the conservation law (3.11). Due to this, and using the dynamic and 
background equations one concludes that zero's component 1° of the conserved current in (3.11), based on 
t a , contains only the first time derivatives of . Therefore 1° has the normal behaviour with respect to 
initial conditions with the definition of the integral quantities analogous to (3.6). Therefore this requirement, 
at least, may be unnecessarily restrictive. 

Babak and Grishchuk give and discuss a wide bibliography of works, where a possibility to consider 
gravitons with non-zero masses is studied [23]. Here, we note only the Visser paper [37]. It is clear that an 
inclusion of non-zero masses of gravitons leads to a non-Einstein gravity. Visser, analyzing foundations of 
GR, find that a more economical way turns out an inclusion of a background metric, and he realizes this 
possibility. Thus a philosophy of the works [23] and [37] coincide. An advantage of the approach in [23] 
is that, using the field formulation of GR, it permits to study the problem from more general positions. 
Thus, the linear Visser's equations are a particular case of Babak and Grishchuk's linear equations, when 
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one sets a = (3. Beginning from the second order, equations in [23] and [37] are different even for a = f3. 
The reason is that Visser defines perturbations as h^ v = <? M „ — g , whereas Babak and Grishchuk use 
in correspondence with our above definitions. It is interesting that the same reason initiates differences in 
definitions of energy-momentum tensors (see section 4). In spite of these differences both of approaches, [23] 
and [37], qualitatively give the same predictions. 
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